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Q ■ Abstract 



Massless matter fields and non-Abelian gauge fields are localized on domain walls 
in a (4+l)-dimensional U(N) C gauge theory with SU(N) L x SU(N) R x U(1) A fla- 
vor symmetry. We also introduce SU(N) L+R flavor gauge fields and a scalar-field- 
dependent gauge coupling, which provides massless non-Abelian gauge fields local- 
ized on the wall. We find a chiral Lagrangian interacting minimally with the non- 
Abelian gauge field together with nonlinear interactions of moduli fields as the (3+1)- 
dimensional effective field theory up to the second order of derivatives. Our result 
provides a step towards a realistic model building of brane-world scenario using topo- 
logical solitons. 
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1. Introduction 



Gauge hierarchy problem is a good guiding principle to construct theories beyond the 
Standard Model (SM). Brane world scenario^ is one of the most attractive proposals to 
solve this problem, besides models with supersymmetry (SUSY).™ In the brane world sce- 
nario, it is assumed that all fields except the graviton field are localized on (3+l)-dimensional 
world volume of a defect called 3-brane, immersed in a many- dimensional space-time called 
bulk. In order to realize such a scenario dynamically, we may use a topological soliton. 
For instance, let us consider a domain wall solution as the simplest soliton. To obtain 
(3+l)-dimensional world volume on the domain wall, we need to consider a theory in a 
(4+l)-dimensional space-time. Bulk fields in (4+l)-dimensions can provide massless modes 
localized on the domain wall, besides many massive modes in general. After integrating 
over massive modes, one obtains low-energy effective field theory describing the effective 
interactions of massless modes. Massless matter fields have been successfully localized on 
domain wallsP but localization of the gauge field on domain walls in field theories has been 
difficult P It has been noted that the broken gauge symmetry in the bulk outside of the soli- 
ton inevitably makes the localized gauge field massive with the mass of the order of inverse 
width of the wallP® To localize a massless gauge field, one needs to have the confining 
phase rather than the Higgs phase in the bulk outside of the soliton. Earlier attempts used 
a tensor multiplet in order to implement Higgs phase in the dual picture, but this approach 
successfully localize only £7(1) gauge fieldP More recently, a classical realization of the con- 
finement!^ 1 ^ through the position-dependent gauge coupling has been successfully applied 
to localize the non-Abelian gauge field on domain walls. 12 ' The nontrivial profile of this 
position-dependent gauge coupling was naturally introduced on the domain wall background 
through a scalar-field-dependent gauge coupling function resulting from a cubic prepotential 
of supersymmetric gauge theories. The appropriate profile of the position-dependent gauge 
coupling was obtained from domain wall solutions using two copies of the simplest model 
or from a model with less fields and a particular mass assignment. However, it was still a 
challenge to introduce matter fields in nontrivial representations of the gauge group of the 
localized gauge field. 

Parameters of soliton solutions are called moduli and can be promoted to fields on the 
world volume of the soliton. Massless fields in the low-energy effective field theory on the 
soliton background are generally given by these moduli fields. Moduli with non-Abelian 
global symmetry is often called the non-Abelian cloud, and has been explicitly realized in 
the case of domain walls using Higgs scalar fields with degenerate masses in U(N) C gauge 
theories.^ This model also has a non-Abelian global symmetry SU (N) L x SU (N) R x U (I) a, 
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which is somewhat similar to the chiral symmetry of QCD. If we turn this global symmetry 
into a local gauge symmetry, we should be able to obtain the usual minimal gauge coupling 
between these moduli fields and the gauge field. Since we wish to localize the gauge field on 
the domain wall, it is essential to choose the global symmetry of moduli fields to be unbroken 
in the vacua (of both left and right bulk outside of the wall). This choice will guarantee that 
the bulk outside of the domain wall is not in the Higgs phase. Therefore we are led to an idea 
where we introduce gauge fields corresponding to a flavor symmetry group of scalar fields 
which will be unbroken in the vacuum. If we introduce the additional scalar-field-dependent 
gauge coupling function similarly to the supersymmetric model, we should be able to localize 
both massless matter fields and the massless gauge field at the same time on the domain 
wall. 

The purpose of this paper is to present a (4+l)-dimensional field theory model of localized 
massless matter fields minimally coupled to the non-Abelian gauge field which is also localized 
on the domain wall with the (3+l)-dimensional world volume. We also derive the low-energy 
effective field theory of these localized matter and gauge fields. To introduce non-Abelian 
flavor symmetry (to be gauged eventually) in the domain wall sector, we replace one of the 
two copies of the U(l) c gauge theory with the flavor symmetry [7(1) £ x U(1)r in Ref JTB"]) . by 
U (N) c gauge theory with the extended flavor (global) symmetry SU (N) L x SU (N)r x U(1)a- 
By choosing the coincident domain wall solution for this domain wall sector, we obtain the 
maximal unbroken non-Abelian flavor symmetry group SU(N)l + r which is preserved in 
both left and right vacua outside of the domain wall. Therefore we can introduce gauge field 
for the (subgroup of) the flavor SU(N)l + r symmetry. In order to obtain the field-dependent 
gauge coupling function, for the gauge field localization mechanism,'® we also introduce a 
coupling between a scalar field and gauge field strengths inspired by supersymmetric gauge 
theories, although we do not make the model fully supersymmetric at present. This scalar- 
field-dependent gauge coupling function gives appropriate profile of position-dependent gauge 
coupling through the background domain wall solution. With this localization mechanism 
for gauge field, we find massless non-Abelian gauge fields localized on the domain wall. 
We also obtain the low-energy effective field theory describing the massless matter fields in 
the non-trivial representation of non-Abelian gauge symmetry. Since our flavor symmetry 
resembles the chiral symmetry of QCD before introducing the gauge fields that are localized, 
we naturally obtain a kind of chiral Lagrangian as the effective field theory on the domain 
wall. We find an explicit form of full nonlinear interactions of moduli fields up to the second 
order of derivatives. Moreover, these moduli fields are found to interact with SU(N)l + r 
flavor gauge fields as adjoint representations. In analyzing the model, we use mostly the 
strong coupling limit for the domain wall sector. The strong coupling is merely to describe 
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our result explicitly at every stage. Even if we do not use the strong coupling, the physical 
features are unchanged. It is easy to expect that (the part of) the gauge symmetry is broken 
when the walls separate in each copy of the domain wall sector. Our results of the low- 
energy effective field theories shows that flavor gauge symmetry SU (N) l+r is broken on the 
non-coincident wall and the associated gauge bosons acquire masses as walls separate. This 
geometrical Higgs mechanism is quite similar to D-brane systems in superstring theory. So 
our domain wall system provides a genuine prototype of field theoretical D3-branes. This is 
an interesting problem, which we plan to analyze more in future. We also find indications 
that additional moduli will appear in the supersymmetric version of our model, which is also 
an interesting future problem to study. 

The organization of the paper is as follows. In section [21 we explain the localization 
mechanism by taking Abelian gauge theory as an illustrative example. In section [31 we 
introduce the chiral model with the non- Abelian flavor symmetry for the domain wall sector 
and then also introduce gauge fields for the unbroken part of the flavor symmetry. By 
introducing the scalar-field-dependent gauge coupling function, we arrive at the localized 
massless gauge field interacting with the massless matter field in a nontrivial representation 
of flavor gauge group. The low-energy effective field theory is also worked out. In section 
HI an attempt is made to make the model supersymmetric. New additional features of 
the supersymmetric models are also described. In section El we summarize our results and 
discuss remaining issues and future directions. In Appendix [A] we discuss domain wall 
solution for gauged massive CP 1 sigma model. Appendix [B] describes derivation of effective 
Lagrangian which includes full nonlinear interactions between moduli fields. Appendix [C] 
contains derivation of positivity condition for the potential appearing in section HJ 



§2. Abelian-Higgs model of gauge field localization 

2.1. The domain wall sector 

Let us illustrate the localization mechanism for the gauge fields and the matter fields on 
the domain walls by using a simplest model in (4+l)-dimensional spacetime : two copies 
{i = 1, 2) of U(l) models, each of which has two flavors (L, R) of charged Higgs scalar fields 
H{ = (H iL , H iR ) : 

C i = "A (Hin) 2 + A ( d M°i? + \VmH1 2 - Vi, (2-1) 

2 

Vi = ^ (\H t \ 2 - v 2 ) 2 + \<TiHi - E t M t \ 2 . (2-2) 
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We use the metric tjmn = diag(+, — , • • • , — ), M,N = 0, 1, • • • ,4. The Higgs field Hi is 
charged with respect to the U(l)i gauge symmetry and the covariant derivative is given by 



V M Hi = d M Hi + iw^Hi, (2-3) 

where w l M is the U(l)i gauge field with the field strength 

Hin = 9mw% - d N wir. (2-4) 

Since we want domain walls, we will choose 

Mi = diag (mi, -mi) , (m» > 0), (2-5) 

resulting in the U(l)iA flavor symmetrjo We have included the neutral scalar fields Oi in this 
Abelian-Higgs model. The gauge coupling appears not only in front of the kinetic terms 
of the gauge fields and <7j, but also as the the quartic coupling constant of Hi. Both these 
features are motivated by the supersymmetry. Indeed, we can embed this bosonic Lagrangian 
into a supersymmetric model with eight supercharges by adding appropriate fermions and 
bosons, which will not play a role to obtain domain wall solutions. We have taken this special 
relation among the coupling constants only to simplify concrete computations below. One 
may repeat the following procedure in models with more generic coupling constants without 
changing essential results. 

The first term of the potential is the wine-bottle type and the Higgs fields develop non- 
zero vacuum expectation values. There are two discrete vacua for each copy i 

(H iL , H iR , Gi) = (vi, 0, mi) , (0, Vi, -m { ) . (2-6) 

Thanks to the special choice of the coupling constants in motivated by the super- 
symmetry, there are Bogomol'nyi-Prasad-Sommerfield (BPS) domain wall solutions in these 
models. Let y be the coordinate of the direction orthogonal to the domain wall and we 
assume all the field depend on only y. Then, as usual, the Hamiltonian can be written as 
follows 

Hi = (<V* + 9- {m 2 - v 2 )) 2 + \V y H t + a t H t - H t M t \ 2 



+ v 2 d y a t - d y ({(TtHi - HM)Hf 
> v 2 d y a t - d y ({(TtHi - HiMi)Hpj . (2-7) 



*' Phase rotation of and Hir in the same direction f/(l), is gauged and the remaining global 
symmetry is in the opposite direction and is denoted as U{\)iA- 
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Thus the Hamiltonian is bounded from below. This bound is called Bogomol'nyi bound, and 
is saturated when the following BPS equations are satisfied 

dy(n + g\ [\H,\ 2 - vf) = 0, V v Hi + aiH, - HM = 0. (2-8) 

In order to obtain the domain wall solution interpolating the two vacua in Eq. (12 -6p . we 
impose the boundary conditions : 

(H iL , H iR , (n) = (0, v u -rrii) , y = -oo, 

(H iL , H iR , (Ji) = (v t , 0, mi) , y = oo. (2-9) 
Tension Tj of the domain wall is given by a topological charge as 

T l= dy vldyVi-dyUaiHi-HiMjHn = 2m^ 2 . (2-10) 

J-oo 1 V /J-oo 

The second equation of the BPS equations f l2-8jl can be solved by the moduli matrix 
formalismpJ'^ with the constant matrix (vector) Hio = (Cn, Cm) 



Hi = Vie-^H i0 e M * y , a { + iw t = ^d y ^. (2-11) 



For a given H i0 , the scalar function ^ is determined by the master equation 

dfa = 2g*v* (l - e-^H^Hl) . (2-12) 

The asymptotic behavior of the field ipi is determined by the condition that the configuration 
reaches the vacuum at left and right infinities: 

if>i^ log H i0 e 2M ^Hj , \y\^oo. (2-13) 

There exists redundancy in the decomposition in Eq. (12-llj) . which is called the ^/-transformation: 

H i0 ^ViH i0 , ^ ^ ^ + 2hgV h ViEC*. (2-14) 

For example, a single domain wall solution centered at y = can be generated by a moduli 
matrix 

H l0 = (l,l). (2-15) 

Then the master equation is 

QjVi = 2 9itf (1 - (e 2m ' y + e- 2m ' y )) . (2-16) 
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Fig. 1. The left panel shows profiles of Hm (solid line), Hm (long-dashed line), and <j{ (dashed- 
line) with finite gauge coupling (gj = 0.5). The right panel shows a plot of a^. dashed curve 
for finite (jji = 0.5) gauge coupling and solid curve for strong gauge coupling (gi = oo). The 
other parameters are mi = Uj = 1. 

No analytic solutions for the master equation have been found for finite gauge couplings gi, 
so we must solve it numerically. The corresponding solution is shown in Fig. [TJ The generic 
solutions of the domain wall are generated by the generic moduli matrices (after fixing the 
^-transformation) 

Hm = (Co,, Cm) , Cn, Cm e C*. (2-17) 

The complex constants Cm, Cm are free parameters containing the moduli parameters of the 
BPS solutions. The moduli parameter can be denned by 

d = = e iai e mm . (2-18) 

The other degree of freedom in Cm, Cm can be eliminated by the V^-transformation in 
Eq. (12-141) and has no physical meaning. Then the master equation is found to be 

d 2 y i)i = 2g 2 v 2 (1 - ( e 2mi{y - m) + e ~ 2mdy - yi) )) . (2-19) 

It is obvious that the real parameter is the translational moduli of the domain wall. 
The other parameter a« is an internal moduli which is the Nambu-Goldstone (NG) mode 
associated with the U(l)iA flavor symmetry spontaneously broken by the domain walls. 

One can take, if one wishes, the strong gauge coupling limit of the Lagrangian Li. As 
is well-known, the U(l) gauge theory with two flavors of Higgs scalars in the strong gauge 
coupling limit becomes a non-linear sigma model whose target space is CP 1 : 

\H t \ 2 = \H lL \ 2 + \H lR \ 2 = vl (2-20) 

The gauge fields and the neutral scalar field become infinitely massive and lose their kinetic 
terms. They are mere Lagrange multipliers in the limit, and are solved as 

< = -^2 (Hid M H\ - d M HiHpj , a t = ^HMHj. (2-21) 
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Plugging these into we get 

C? = d M HiPid M H} - HtMiPMHl (2-22) 

with a projection operator 



1 



Pi = 1 - ^HjHi. (2-23) 



Let us introduce an inhomogeneous coordinate fa of CP 1 by 



H lL = - Vi = , H iR = - V jt =. (2-24) 
VTTW y/T+W 

Then the Lagrangian of the CP 1 model in terms of fa is 

cr = v fMfa?-^\fa\^ 

(i + \fa\ 2 ) 2 

Let us reconsider the domain wall solutions in this limit. The Hamiltonian can be written 

as 

n? = (i + ijfeiV + 2m ' '' 2 + 2m,t ' TyTTW 



>- 2m -"-^TTW- (2 ' 26) 

The BPS equation and the boundary conditions are given by 

dyfa + 2mifa = 0, 

fay = -oo) = oo, fay = oo) = 0, (2-27) 

corresponding to the boundary conditions in Eq. (12-9|) . The BPS equation can be easily 
solved by 

fa = C^C iR e- 2miV = C?e- 2m > y . (2-28) 
The tension of the domain wall is 

POD 

dy Uf = 2m iV 2 . (2-29) 



This is the same as the one in the finite gauge coupling model. 

In this way, the strong gauge coupling limit has a great advantage compared to the finite 
gauge coupling case. One can exactly solve the BPS equation and see the moduli parameter 
in the analytic solutions. Furthermore, there is no important differences between domain wall 



8 



solutions in the finite coupling (Abelian-Higgs model) and the strong coupling (non-linear 
sigma model). Both solutions have the same tension of domain wall and the same number 
of the moduli parameters. To see the difference explicitly, let us compare the configuration 
of the neutral scalar field 0{. In the strong gauge coupling limit, it can be written as 



<Ti=m } : rn? = mi tanh2mi(y - yi), (2-30) 



where we have used 



C i = e ia e miVi . (2-31) 

In Fig. (U we show the configurations of o\ in two cases, the one in the small finite gauge 
coupling and the one in the strong gauge coupling limit. As can be seen from the figure, 
there are no significant differences. 

Let us next derive the low energy effective theory on the domain wall. We integrate all the 
massive modes while keeping the massless modes. We use the so-called moduli approximation 
where the dependence on (3+l)-dimensional spacetime coordinates comes into the effective 
Lagrangian only through the moduli fields: 

d -> Ci(x»), Uv) -> <f>i{y,Ci(x»)) = c^fe- 2 ^. (2-32) 

The effective Lagrangian for the moduli field Cj(x M ) can be obtained by plugging this into 
the Lagrangian Li and integrate it over y. This can be done explicitly as follows. 
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Aeff= / dy ^ rZU = ? 9 ■ I 2 " 33 ) 



With Eq. (12-311) . the effective Lagrangian is given by 

k# = ^(d,y t f + ^{d, ai )\ (2-34) 

where energy of soliton solution is neglected since it does not contribute to dynamics of mod- 
uli. Note that 2rriivf is precisely the domain wall tension. This is the free field Lagrangian. 

Although we have derived this effective Lagrangian in the strong gauge coupling limit, 
we can obtain the same Lagrangian in the finite gauge coupling constant. In other words, 
the effective Lagrangian cannot distinguish the infinite versus finite coupling cases at least 
in the quadratic order of the derivative expansion. 
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2.2. Localization of the Abelian gauge fields 

In the previous subsection, we have seen that the NG modes of the translation and U(l) 
global symmetry are the only massless modes in the Abelian-Higgs model. They are localized 
on the domain wall. There are no massless gauge field on the domain wall and all the modes 
contained in the gauge field are massive. The mass of the lightest mode of the gauge field 
is of the order of the inverse of the width of the domain wall, since the bulk outside of the 
domain wall is in the Higgs phase. The low energy effective Lagrangian for the massless 
fields is obtained after integrating out the massive modes including gauge fields. 

In order to obtain the massless gauge field to be localized on the domain wall, we need 
a new gauge symmetry which is unbroken in the bulk. Recently, a new mechanism was 
proposed to localize gauge fields on domain walls 

A key ingredient is the so-called dielectric coupling const ant?® ^ for the new gauge 
symmetry. To illustrate the new localization mechanism, let us introduce a new U(l) gauge 
field <2m which we wish to localize on the domain wall. Since this gauge symmetry should 
be unbroken in the bulk, we consider the case where all the Higgs fields are neutral under 
this newly introduced U(l) gauge symmetry. The gauge field au is assumed to couple to 
the neutral scalar fields oi only in the following particular combination 

C = & + £2 - $ ( — - —) (Gain) 2 , (2-35) 
2 \mi 1712/ 

where a real constant A with the unit mass dimension, in accordance with the (4+1)- 
dimensional spacetime and the field strength is defined by 

Gain = Qm^n — Qn^m- (2-36) 

The field-dependent gauge coupling function is given by 

1 _ A / ox o~ 2 
4e 2 (cr) 2 \m\ m2 

which depends on the position y through fields (7j. Thus the field-dependent gauge coupling 
function e(a) plays the role of the dielectric coupling constant. Furthermore, the special 
choice in Eq. ( I2-37P is chosen for the gauge interaction to become strongly coupled in the 
bulk (crj — > ±rrii as y — > ±00). 

Let us again consider a double copy of domain walls as a background configuration in the 
Abelian-Higgs model in Eq. (!2-35[) . Since Lagrangian has no term linear in om, the equations 
of motion for clm is trivially solved by om = 0, and the rest of the equations of motion are 
explicitly the same as those in the previous subsection. Therefore the domain wall solution 



(2-37) 
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in the previous subsection together with = is still a solution of the equations of motion. 
Clearly, the low energy effective Lagrangian on the domain wall is also unchanged 



£ c ff — A.eff + ^2,eff — -^(GfiuY 



(2-38) 



except for the additional kinetic term (the last term) of the (3+l)-dimensional gauge field 
Wfj,, which is the zero mode (y-independent mode) of the (4+l)-dimensional field w^. The 
(3+l)-dimensional gauge coupling constant is given by 



1 A 

<±e 4 z, >/ _ 00 



dy 



&2 \ 


A 




^2 


m 2 ) 


" 4 


mi 


m 2 



Kv2 - yi) 



(2-39) 



where we have used the asymptotic behavior ipi — > log 2 cosh 2mj (y — y.\) as \y\ — )■ oo. Note 
that this result is again independent of the gauge couplings gi in the domain wall sector. In 
summary, the low energy effective Lagrangian is 

,2 „,2 



c 



cff 



8=1,2 



—z—idnyi) +—(drati) 

Z v/7,,; 



- A(y 2 - yi)(C? ( 



A 1 ^ y 



(2-40) 



Now we separate the quantum fields (fluctuations) from the classical background moduli 
parameters by 



yi{x^) = y® + 8y u a^) = a? + 5a { . 



(2-41) 



Then the effective Lagrangian up to the second order of the small quantum fluctuations is 
given by 



i=l,2 



2m,t> 2 



rrii 



Ky°2 



vl){G. 



[IV j 



(2-42) 



We note that the massless gauge field has a positive finite gauge coupling squarec0 
l/(4A(y 2 ° - ?/?)) provided y° 2 - y\ > 0. 

Although we succeeded in localizing the massless U(l) gauge field on the domain 
walls, the Lagrangian Eq. f)2-42p has no charged matter fields minimally coupled with the 
localized gauge field a /t . To obtain matter fields interacting with the localized gauge field, 
one may be tempted to identify the Higgs fields Hi = (H iL , H iR ) as matter fields 



**) 



with 



*) Here we are content with the fact that the positivity of the gauge kinetic term is assured at least in 
finite region of moduli space, instead of just at a point. However, it is possible to make a more economical 
model where one has less moduli, and the positivity of the gauge kinetic term is assured!^ 1 

**> We consider the diagonal subgroup U(1)a of U(1)ia and U(V)2A ■ Actually the U{\)iA global sym- 
metries are broken by the domain wall solution, we consider this gauging to leading order of gauge coupling 
only to illustrate the Higgs mechanism for the broken symmetry. 
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charges (1,-1). The minimal gauge interaction of Higgs fields with the au is introduced 
through the modified covariant derivatives as 

i> M Hi L = d M H lL + iw l M H iL + ia M H iL , (2-43) 
VmHir = d M H iR + iw l M H iR - ia M H iR . (2-44) 

Since the moduli field Cj is charged, the derivatives in the low energy effective theory 
Eq. (I2-33P should be replaced by the covariant derivative 

d^d ->• VyPi = d^a + ia^d. (2-45) 

It is straightforward task to derive the effective Lagrangian with the covariant derivative 
above along the same line of reasoning for the previous case 

£e ff (% ,a°)=E 

i=l,2 

-Ky 2-y°i)(Q^ 2 - (2-46) 

This clearly shows that the new gauge field a M is not massless due to the Higgs mechanism, 
and should be integrated out together with the other massive fields. Namely the low energy 
effective Lagrangian does not include the massless gauge fields, since the U(l) symmetry 
which we gauged is broken by the domain wall. A more explicit example at the strong gauge 
coupling limit is described in Appendix |A] Thus the Abelian-Higgs model in this section 
gives an important lesson that we should not gauge a symmetry which is broken by the 
domain wall solution, since the corresponding gauge fields may be localized on the domain 
walls but they become massive and should be integrated out from the low energy effective 
theory. In the next section, we will give a model with a non-Abelian global symmetry whose 
unbroken subgroup can be gauged to yield massless localized gauge fields on the domain 
wall. 



Uli 



(djai + qia h 



§3. The chiral model 

In this section we study domain walls in the chiral model which is a natural extension of 
the Abelian-Higgs model in the previous section. This chiral model leads to two important 
consequences 1) massless non-Abelian gauge fields are localized on the domain wall and 
moreover 2) the scalar fields which are non-trivially interacting are also localized on the 
domain walls. 



12 



3.1. The domain walls in the chiral model 

As a natural extension of the domain wall sector in the previous section, we consider the 
Yang-Mills-Higgs model with SU(N) C x U(l) gauge symmetry with S[U(N) L x U(N) R ] = 
SU(N)l x SU(N)r x U(1)a flavor symmetry.^ JS) To localize the gauge field in a simple 
manner, we again introduce two sectors £i and £2, but only the former is extended to Yang- 
Mills-Higgs system and the latter is the same form as (12 -ip . The second sector couples to the 
first sector through the coupling as described in (I2-35P after gauging the flavor symmetry 
it plays a role as localization of gauge fields, combined with the first sector. The matter 
contents are summarized in Table [B Since the presence of two factors of SU(N) global 
symmetry resembles the chiral symmetry of QCD, we call this Yang-Mils-Higgs system as 
the chiral model. 





SU(N) C 


U(l)i 


U(l)2 


SU{N) L 


SU{N) R 


U(l) 1A 


U(1) 2A 


mass 
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1 





□ 


1 


1 





m x l N 


HlR 


□ 


1 





1 


□ 


-1 





-miljv 


Hi 


adj © 1 








1 


1 











H2L 


1 





1 


1 


1 





1 


m 2 


H2R 


1 





1 


1 


1 





-1 


-m 2 


s 2 


1 








1 


1 












Table I. Quantum numbers of the domain wall sectors in the chiral model. 
The Lagrangian is then given by 



1 



£ = d + £ 2 
£1 = Tr 

Vi = Tr 



2gi 



1 

7i 



^^IMiv) 2 H 2^ M ^^) 2 + \PM~Rl\ 



V 1} 



(3-1) 
(3-2) 

(3-3) 



with Hi = (Hil, H 2 l)- £2 is the same form as (12- lj) with i = 2. Gauge fields of U(N) C = 
(SU (N) c x U (1)i)/Zn are denoted as Wim, and adjoint scalar as E\. The covariant derivative 
and the field strength are denoted as V M Ei = 8mEi + i \Wim > ^1] > T^mH\ = 3mHi + 
iWi M Hi, and F 1MN = d M Wi N — d N Wi M + i [Wim, Win]. The mass matrix is given by 
Mi = diag (miljv, — miljv). Let us note that the chiral model reduces to the Abelian-Higgs 
model in the limit of iV — » 1, by deleting all the SU(N) groups. 

The second sector is just necessary to realize the field- dependent gauge coupling function 
similar to f !2-35|) as we will discuss in the subsequent subsection. In the rest of this subsection, 
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we focus only on the first sector (i — 1) and suppress the index % — 1. The symmetry 
transformations act on the fields as 



H=(H L ,H R )^U C (H L ,H R ) 
E -> C/^C/], 



U R e 



(34) 
(3-5) 



with C/ c G C/(7V) C , l/ L G 5C/(7V) L , U(N)r G 5C/(7V) r and G 

There exist A^ + 1 vacua in which the fields develop the following VEV 



H=(H L ,H R )=v 

IjV-r 



Ljv-r 



0, 







N-r 



m 



(3-6) 
(3-7) 



with r = 0, 1, 2, • • • , N. We refer these vacua with the label r. In the r-th vacuum, both 
the local gauge symmetry U(N) C and the global symmetry are broken, but a diagonal global 
symmetries are unbroken (color-flavor-locking) 



U(N) C x SU(N) L x SU(N) R x U(1) A -> 

ST/(jV - r) L+c x SC/(r) L x SU(r) R+c x 5C/(7V - r) R x C/(1) A+C 



(3-8) 



As in the Abelian-Higgs model, the BPS equations for the domain walls can be obtained 
through the Bogomol'nyi completion of the energy density with the assumption that all the 
fields depend on only the fifth coordinate y and = 0: 

2 



H = Tr 



1 ( 2 \ 

-(v y E- 9 -( V 2 l N -HH^ + \V y H + EH — HM\ 



+ d y {Tr [v 2 E - (EH - HM) ff f ] } 
> d y {Tr [v 2 E - (EH - HM) ff f ] } . 

This bound is saturated when the following BPS equations are satisfied 

V y E - 9 - (v 2 l N - HH^) = 0, 
V y H + EH- HM = 0. 

The tension of the domain wall is given by 

/oo 
dy d y {Tr [v 2 E - (EH - HM) if t] } 
-oo 

= v 2 Tr [27(+oo) - E(-oo)} . 



(3-9) 



(3-10) 
(3-11) 



(3-12) 
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Let us concentrate on the domain wall which connects the 0-th vacuum at y — > oo and 
the iV-th vacuum at y — > — oo. Its tension can be read as 



T = 2Nv 2 m, (3-13) 

from Eq. (13- 12p . Since there are N + 1 possible vacua, the maximal number of walls is 
TV at various positions. The simplest domain wall solution corresponding to the coincident 
walls is given by making an ansatz that H^, H R , E and W y are all proportional to the unit 
matrix. Then the BPS equations (I3-I0p and (13-lip can be identified with the BPS equations 
in Eq. (12 • 8 1) in the Abelian-Higgs model. Thus the domain wall solution can be solved as 

H L = ve-ie my 1 N , (3-14) 
H R = ve-*e- mv 1 N , (3-15) 

Z + iW v = ~d v il>l N , (3-16) 

where ip is the solution of the master equation (12- 12|) in the Abelian-Higgs model. Eq. ()3-8p 
shows that the unbroken global symmetry for iV-th vacuum {Hi = 0, H R = vljy and 
E = — mix) at the left infinity y —> — oo is SU(N)l x SU(N)r +c x U(1)a+c, whereas that 
for the 0-th vacuum (Hl = vIn, Hr = and E = mix) at the right infinity y — » oo is 
SU(N) L+C x SU(N) R x U(1) A+C . 

The domain wall solution further breaks these unbroken symmetries because it interpo- 
lates the two vacua. The breaking pattern by the domain wall is *3 

U(N) C x SU(N) L x SU(N) R x U(1) A SU(N) L+R+C . (3-17) 

This spontaneous breaking of the global symmetry gives NG modes on the domain wall as 
massless degrees of freedom valued on the coset similarly to the chiral symmetry breaking 
in QCD : 

SUm, x SUiN)* x 



SU(N) 



L+R+c 



Since our model can be embedded into a supersymmetric field theory, these NG modes (U(N) 
chiral fields) appear as complex scalar fields accompanied with additional iV 2 pseudo-NG 
modes **). 



*) The unbroken generators of U(1)a+ c for r-th vacuum contains different combination of U(N) C gen- 
erators depending on r. Therefore the right and left vacua preserve actually different U(1)a+c, and the wall 
solution does not preserve any of these U(1)a+c- 

**) One of them is actually a genuine NG mode corresponding to the broken translation. 
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3.2. Localization of the matter fields 

In the remainder of this subsection, we will give the low-energy effective Lagrangian on the 
domain walls where the massless moduli fields (the matter fields) are localized. The best way 
to parametrize these massless moduli fields is to use the moduli matrix formalism^ 1 ^ 

H L = ve my S- 1 , (3-19) 
H R = ve^S^e* , (3-20) 
Z + iWy = S^dyS, (3-21) 

where S G GL(N, C) and ft = SS^ is the solution of the following master equation 

d y (ft^dyft) = g 2 v 2 (1 N - ft- l ft ) , (3-22) 

where 

ft = e 2my l N + e - 2m Ve^ . 

We have used the ^-transformation to identify the moduli e^, which is a complex by 
matrix. It can be parametrized by an iV x iV hermitian matrix x and a unitary matrix U 

= e^\j\ (3-23) 

where U is nothing but the U(N) chiral fields associated with the spontaneous symmetry 
breaking Eq. (13-18P and x is the pseudo-NG modes whose existence we promised above. 

In the strong gauge coupling limit g — > oo, solution of master equation is simply ft = ft . 
After fixing the U(N) C gauge, we obtain 

S = e £/ V2 cosh(2rm/ - x) . (3-24) 

Let us denote, for brevity 



y = 2my — x , (3-25) 



the Higgs fields are then given as 



e y/2 



H L = v^== , (3-26) 

V 2 cosh y 

p-fi/2 

H R = v -==rf . (3-27) 

V 2 cosh y 

From this solution, one can easily recognize that eigenvalues of x correspond to the positions 

of the Af domain walls in the y direction. Now we promote moduli parameters x and U to 
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fields on the domain wall world volume, namely functions of world volume coordinates x M . 
We plug the domain wall solutions HL jR (y;x(x tJ, ),U(x ,x )) into the original Lagrangian C in 
Eq. (13-21) at g — > oo and pick up the terms quadratic in the derivatives. Thus the low energy 
effective Lagrangian is given by 

/oo 

dyTi d il H L d^H\ j + d il H R d l "H R -v 2 W lx W 11 , (3-28) 
-oo 

where 

d p H L H[ - H L d„Hl + (L « R)] . (3-29) 



w - = h 



Here we have eliminated the massive gauge field by using the equation of motion. Using 
the solutions for Hi and H R we have found a closed formula for the effective Lagrangian up 
to the second order of derivatives but with full nonlinear interactions involving moduli fields 
x and U. Detailed derivation is given in Appendix [Bj 

Here we exhibit the result only in the leading orders of U — 1 and x: 

C cS = — Tr (djfldPU + d.xd^x) + ... (3-30) 
2m V / 

When N = 1 and with the redefinitions U = e 2iai , and x = 2my 1 , this coincides with the 

effective Lagrangian £i=i >e ff in Eq. (12-341) of the Abelian-Higgs model, which we obtained in 

the previous section. 

3.3. Localization of the gauge fields 

Let us next introduce the gauge fields which are to be massless and localized on the 
domain walls. As we learned in section El the associated gauge symmetry should not be 
broken by the domain walls. Therefore, the symmetry which we can gauge is the unbroken 
symmetry SU(N) L+R+C itself or its subgroup. 

Let us gauge SU(N) L+R = SU(N) V and let be the SU(N) L+R gauge field. The Higgs 
fields are in the bi-fundamental representation of U(N) C and SU(N)l +r . The covariant 
derivatives of the Higgs fields are modified by 

V M H 1L = d M H lL + iW 1M H 1L - iH 1L A M , (3-31) 
V M H lR = d M H 1R + iW lM H 1R - iH lR A M . (3-32) 

The quantum numbers are summarized in Table [Til 

We now introduce a field-dependent gauge coupling function g 2 (S) for Am, which is 
inspired by the supersymmetric model in Ref JT2l) . 

_L_ = * (T^l _ ^l) (3 .33) 
2e 2 (r) 2 \Nmi m 2 J ' 1 } 
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bU(N) c 




T T / 1 \ 


OT T f A 7"\ 


T T / 1 \ 


U{l) 2 A 


mass 




□ 


1 





□ 


1 





mil N 


tl\R 


1 — 1 

U 


1 





1 — 1 

u 


— 1 





-mil^ 


Zi 


adj © 1 








1 













1 





1 


1 





1 


m 2 




1 





1 


1 





-1 


— m 2 


Z2 


1 








1 












Table II. Quantum numbers of the domain wall sectors in gauged chiral model 



The Lagrangian is given by 

£ = A + C 2 - ^^yTr [G MN G MN ] . (3-34) 

The Ci in Eq. f)3-34p is given by Eq. f)3-2p where the covariant derivatives are replaced with 
those in Eqs. (ESI]) and (13^32]) . 

We first wish to find the domain wall solutions in this extended model. As before, we 
make ansatz that all the fields depend on only y and = = 0. Let us first look at the 
equation of motion of the new gauge field Am- It is of the form 

V M G MN = J N , (3-35) 

where Jm stands for the current of A M . Note that the current J M is zero, by definition, if 
we plug the domain wall solutions in the chiral model before gauging the SU(N)l+r. This 
is because the domain wall configurations do not break SU(N)l + r. Therefore, Am = is a 
solution of Eq. f)3-35p . 

Then, we are left with equation of motion with Am = which are identical to those in 
the ungauged chiral model in the previous subsections. Therefore the gauged chiral model 
admits the same domain wall solutions as those (Eqs. ( I3-26P and (I3-27P ) in the ungauged 
chiral model. 

The next step is to derive the low energy effective theory on the domain wall world- 
volume in the moduli approximation as in the previous subsections. Again, we promote the 
moduli parameters as the fields on the domain wall world-volume and pick up the terms up 
to the quadratic order of the derivative d^. Similarly to section 13. 2[ we utilize the strong 
gauge coupling limit gi — > 00, to simplify the computation without changing the final result. 
Let us emphasize that we keep the field-dependent gauge coupling function e{E) finite. 
The spectrum of massless NG modes is unchanged by switching on the SU(N) L+R gauge 
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interactions 



We just repeat the similar computation to those in section [3721 Again we shall focus on 
the first sector C\ and suppress the index i = 1 of fields. Since color gauge fields becomes 
auxiliary fields and eliminated through their equations of motion, it is convenient to define 
the covariant derivative only for the flavor (SU (N) l+r) gauge interactions as 



(3-36) 



Then we obtain the effective Lagrangian of the first sector as 



efi 



dy Tr D^H L (D»H L y + D^H R (D»H R y - v'W^W 



2e 2 (E) 



with 



2v 2 



b^H L H{ - H L (D^H L y + (L^R) 



(3-37) 



(3-38) 



Eliminating W^, we obtain the following expression for the integrand of the effective La- 
grangian after some simplification 



-y poo 



2v< 



dy Tr 



(3-39) 



where we defined fields H a b with the label ah of adjoint representation of the flavor gauge 
group SU (N) l+r+c and the covariant derivative as 



T^^H ab = d^Hab + i[Afj,, Hob], H ab = H^Hb, a,b = L, R. 



(340) 



In Appendix [Bj we will describe fully the procedure to derive the effective Lagrangian by 
substituting (13-261) and (13-271) and rewriting in terms of moduli fields x and U. Here we 
merely state the result: 

2 r cosh(£ £ ) — 1 f 1 + tanh^ 



C-[ off — — Tr 

' 2m 



1 — tanh(£. 
t cosh(A) - 1 j ( l + tanh(£,) \ 
+ U U " U £,sinh(£ £ ) ^U-tanKA)/^ X) 

2 tanh(£ £ ) \1 - tanh(£ 4 )y v ; 



(T>^x) 



(341) 



*) Tree level mass spectra are unchanged even though the chiral symmetry SU(N)l x SU(N)r is broken 
by the SU(N)l+r gauge interactions. 
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where 

C A (B) = [A,B] (342) 
is a Lie derivative with respect to A. The covariant derivative is defined by 

V,JJ = dJJ + i[A^U]. (343) 

The above result suggests that the chiral fields U (x M ) and hermitian fields x{x^) are in the 
adjoint representation of SU(N)l + r. Let us now examine the transformation property of U 
and x under the SU(N)l+r flavor gauge transformation on the domain wall background in 
order to demonstrate that they are in the adjoint representation. The domain wall solution 
only preserves the diagonal subgroup SU(N) L+R+C . Eqs.f l34l) and ( 13 • 5[) shows the fields 
transform under the SU (N) l+r+c transformations U as 

H' L = UH L U\ E' R = UE R U\ E' = UEU ] . (344) 

Eqs.([349D and f l3~20|) show that 

S' = USU\ = Ue*U\ Q' = UQU ] . (345) 

The complex moduli is decomposed into hermitian part e x and unitary part U in Eq. fl3-23p . 
Since we can express e 2x = e^e^\ and U = e~^e x , we find that they transform as adjoint 
representations 

By expanding (I34ip . we here illustrate nonlinear interactions of x up to fourth orders in 
the fluctuations x and U — 1 



2 

£ij = ^-TrfvjjtWU + V.xVx + rfVJJ [x, V*. 
2m V 



-i [V„x, x] [x, V»x\ + hpjflU, x\ [x, rfWU] + •••). (347) 

Similarly to Eq.( l2-39|) . we can define the (3+l)-dimensional non-Abelian gauge coupling 
e± by integrating (I3-33P and find 



where yi is the wall position for the z-th domain wall sector. Summarizing, we obtain the 
following effective Lagrangian 



£ e ff — ^i,eff + £>i ,efr — 77-0 Tt" 



(349) 
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where C2, e s is given in ( !2-34[) . This is the main result of this paper. We have succeeded 
in constructing the low energy effective theory in which the matter fields (the chiral fields) 
and the non-Abelian gauge fields are localized with the non-trivial interaction. We show the 
profile of "wave functions" of localized massless gauge field and massless matter fields as 
functions of the coordinate y of the extra dimension in Fig. [2] 
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Fig. 2. The wave functions of the zero modes. DW1 and DW2 stand for the wave functions of the 
massless matter fields of the i = 1 domain wall and i = 2 domain wall, respectively for strong 
gauge coupling limit gi = oo and m» = 1. The gauge fields are localized between the domain 
walls. 

As is seen from Eq. (!3-47[ ) . the flavor gauge symmetry SU(N)l + r +c is further (partly) 
broken and the corresponding gauge field becomes massive, when the fluctuation cj> = e x U 
develops non-zero vacuum expectation values. Especially, x is interesting because its non- 
vanishing (diagonal) values of the fluctuation has the physical meaning as the separation 
between walls away from the coincident case. For instance, if all the walls are separated, 
SU(N) L+R+C is spontaneously broken to the maximal U(l) subgroup {/(l)-^ -1 . However, 
if r walls are still coincident and all other walls are separated, we have an unbroken gauge 
symmetry SU(r) x U(l) N ~ r+l . Then, a part of the pseudo-NG modes x turn to NG modes 
associated with the further symmetry breaking SU(N) L+R+C — >■ SU(r) x U(l) N ~ r+1 , so that 
the total number of zero modes is preservecP^0. These new NG modes, called the non- 
Abelian cloud, spread between the separated domain wallsP^l 1 The flavor gauge fields eat 
the non-Abelian cloud and get masses which are proportional to the separation of the domain 
walls. This is the Higgs mechanism in our model. This geometrical understanding of the 
Higgs mechanism is quite similar to D-brane systems in superstring theory. So our domain 
wall system provides a genuine prototype of field theoretical D3-branes. 



*' In ReffTSf the authors argued that the non-Abelian clouds spreading between walls become massive 
contrary to the results of Refill?]). 
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§4. Embedding into supersymmetric theory 



A crucial point to localize gauge field around domain wall is the coupling between scalar 
and gauge kinetic term. Such a coupling is naturally realized in (4+l)-dimensional super- 
symmetric gauge theory™ This theory generally consists of hypermultiplet part and vector 
multiplet part. The latter is specified by the so-called prepotential. In (4 + l)-dimensional 



serves interactions among vector mutiplets such as fl3-33j) . 

4.1. Supersymmetric model 

In embedding the model into supersymmetric gauge theories in (4 + 1) dimensions, we will 
give non-Abelian global flavor symmetry SU(Ni)v for each copy (i = 1, 2) of the domain wall 
sector, instead of only one copy as in f)3-34p of the previous section. This contains the model 
(l3-34p as a limiting case of N 2 — > 1, and may offer more general situation phenomenologically. 
To formulate supersymmetric gauge theories, we need to introduce Y{ as auxiliary fields 
of U(Ni) c vector multiplet, and <Pj and 3^ as adjoint scalar fields and auxiliary fields of 
SU(Ni) v vector multiplet. As bosonic fields of theories with eight supercharges, we also 
need to double the scalar fields Hi, by introducing another set Hj = (Hj L , H] r ) with masses 
(mjlTVi, — mjlTvJ- They are in the same representations as Hi under U(Ni) c and U{l)iA- 
Explicit charge assignments for hypermultiplets matter fields and adjoint scalar fields are 
summarized in Table IIHI The resultant supersymmetric Lagrangian is written as 



theory the prepotential generally allows up to cubic terms in vector multipletsp^ which 



U(Ni) c U(l) iA SU(Ni) v mass 



HiL 

HiR 

HiL 

Si 



□ ID rriil N . 
DID -niil N% 

□ -1 □ TOjljVi 

□ -1 □ -rriil Ni 
adj 1 



1 adj 



Table III. Quantum numbers of hypermultiplets (Hi, Hi), Si and 




a 



2 



+ Tr { (v M H lL V M Hl L + V M H lL V M H\ L + (L «-> #)) 



i=l 




(4-1) 



22 



where 



C lY = Tr [Hj L Y t H tL - E\ L E lL y % - H^Hl + y i H lL H\ L + (L R) 
V iF = Tr [\EiHiL - H iL {<f> % + mi l Ni )\ 2 + \H iL Si - ($< + mi l N jH iL \ 2 
+\EiH iR - H iR (<Pi - mi l Ni )\ 2 + \H iR Si - & - m l l Nt )H iR \ 2 



(4-2) 



(4-3) 



where a, (3 ■ ■ ■ denote all gauge groups and their generators collectively. We label them with 
the ordering 



a, 13 = O l ,I 1 ,A 1 ;0 2 ,h,A 2 



(4-4) 



where 0, denotes U(l)i parts of U(Ni) c gauge group, while Jj = 1, ■ ■ • ,N 2 — 1 are color 
indices of SU(Ni) c and = 1, • • • , N 2 — 1 denotes flavor indices of SU(Ni) v gauge group. 
The scalar fields S a and auxiliary fields Y a are explicitly written by 



S a = (S Ql , S h , <P Al ; S° 2 , S h ,<P A ' 2 ), 
Y a = (Y 0l , Y h , y Al ; Y 02 , Y h , y M ) , 

and similarly the field strength F^ N and gauge field W" 4 are written by 



MN 



ZJiOl 771/1 /^yAi . 771O2 771/2 \ 

r A/Af) MNi MNi MNi MNi ^MNJ 

Wm = (W° M \W£, A%; <h < 2 ) • 
We adopt the convention of U(Ni) c and SU(Ni)y matrices such as 



Si = S 0i l N . + S h T h , Tr(T 7i T Ji 

1 



<5 = $ Ai T A 



AiiriBi 



Tr(T *T 



S AiB \ (no sum for i) 



(4-5) 
(4-6) 

(4-7) 
(4-8) 



(4-9) 
(4-10) 



Covariant derivatives for H iF and H iR are given as (13-311) and (13 -321) with identical definition 
for Hj L H iR - Covariant derivatives of S 1 * ,<£> Ai are defined as the adjoint representation. We 
will not display the Chern-Simons term Cics and the fermionic term £;f er mion, since we do 
not need them for our analysis. 

Functions a a p{S) are gauge coupling functions, which are given as second derivative of 
the prepotential 

2 



m = E 



i=l 



2W (r ° i)2 + i ( ^ )2 + - 



Xj fs° 

2 



rrii 



JJ0 2 

m 2 



a>ap{S) 



d 2 a(S) 
dE a dEP ' 



(4-11) 



(4-12) 
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From the above prepotential, we see the coupling constants of U(l)i and SU(Ni) c are given 
by §i and gi, respectivel}0. We denote the coupling function of SU(N,i) v corresponding to 
E a = and = $ Bi as e^S), 



em 



mi rri2 



(4-13) 



but will suppress the argument £ to write in the following. 

The constants c a are coefficients of the Fayet-Iliopoulos (FI) terms, allowed to be non-zero 



only for the U(l) part of the gauge groups**) 



Y Y a = c 0l Y°> + c 02 Y° 2 ■ 



(4-14) 



We have assumed both the FI parameters c 0l and c 02 to be positive in the same direction in 
SU(2)r, which is chosen to be along the third component. In this setup, the H fields will 
vanish in the classical solution. Moreover, they do not contribute to the desired order of 
effective Lagrangian. Similarly we have neglected the auxiliary fields Y other than the third 
component in SU(2) R , and we have denoted as Y a . Hence we can call the potential after 
eliminating the auxiliary fields K's to be D-term potential. 

The F-term potential V^p can be worked out from the following superpotential 



Wi = Tr 



{SiHiL - H iL (<P + mi)}H iL + {SiH lR - H iR {$ - m*)} H iR 



(4-15) 



where we restored the tilde fields if's to facilitate writing the superpotential. After elimi- 
nating the auxiliary fields F's, and with the use of 



ViF — —C>iF — WiL\ 2 + \FiR\ 2 + \FiA 2 + |-F1r| 2 . 



(4-16) 



we have (14- 3p . 

Finally, let us work out explicit forms of the D-term potential Vp. Collecting terms 
containing the auxiliary fields F's, we obtain 



i=i 



tL (Y°>) 2 + ^{Y h f + r^Y Ji + (r°< - c 0i )Y°> + ^(y^f + s A *y A > 



+a, ( — - — $ A *y M 

771 1 777,2 



(4-17) 



*> The U(l)i coupling is in principle unrelated to the SU(Ni) c coupling. In section [31 we made a 
simplifying assumption gi = g^f \/2Nl, which allows simple solutions. 

**) The FI parameters CQ i are related to the parameters vf in Eqs. (l3Tp - (|3-3[) in section [3] as c 0i = Nivf . 
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where 



n = H iL H ] iL - H- L HiL + H iR H ] iR - H\ R H iR 



s i = -H] L H iL + H iL H iL - H iR H iR - H iR H iR , 
are Hermitian matrices, with the decomposition 



(4-18) 
(4-19) 



1 

Ni 



(4-20) 



We observe, that in the potential (I4-17P , Y Ti do not couple to the rest of auxiliary fields and 
can be easily eliminated. Having this done, we collect the U(l)i and SU(Ni)y terms into a 
matrix form labeled by a, (3 = 0i, O2, A±, A 2 



Vd = -\j^9Kr U f + \G afs Y«Y? + (r - c) a Y< 



i=l 



(r - c) 0i = r°> - c 0t , (r - c) Ax = s Al . 
Eliminating remaining auxiliary fields we obtain: 



Matrix G = (G a p) is explicitly given by 



G 



( i 





mi 


A2 gvAa 

mi 





1 

If 


m 2 


_ A2_^A 2 

m 2 


mi 


m 2 







\mi 


m 2 








with the inverse 



G- 1 



— x 



/#? - glgWi® 2 -glg 2 m im2 i> 2 

-g\g 2 2 m x m 2 ^ 2 g 2 - g 2 gjm 2 $ 2 



-g{e\m 2 $ Bx 



1 - g 2 <P 2 

-g 2 eim 2 ^ Al 
g\e x m\$ Al 
g 2 2 e x m x $ Bx e\5 MBl - e\g 2 $ AlBl 



\ -g\e 2 m 2 <$ B2 g 2 e 2 mi^> B2 g 2 eie 2 <P Al <P B2 



-g 2 e 2 m 2 <P A2 \ 
g\e 2 m\& A2 
g 2 eie 2 $ M $ Bl 

e | 5 A 2 B 2 _ e 2g2$A 2 B 2 J 



(4-21) 



(4-22) 



(4-23) 



(4-24) 



(4-25) 
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where we abbreviated: 



g 2 = gim 2 2 + g 2 2 mi, (4-26) 
mim2 

$ 2 = $ A i<£ A i + ^ A 2^A 2 ^ (4-28) 

(pAiBi _ (jpfiAiBi _ (pAi^Bi _ (4-29) 



4.2. Positivity of Potential 

The F-term potential (14-30 is manifestly positive. The D-term potential ( 14-210 is positive 
definite under certain conditions. To find the condition we shall decompose ( 14-210 to: 

V D = V 1D + V 2D , (4-30) 
V 1D = \gl{r I ^ + \gl{r 1 >)\ (4-31) 
V 2D = -(G^^r - c) a (r - c) p . (4-32) 

It is clear that the Vm is positive definite by itself. Therefore we can only focus on V2D, 
which is positive if and only if G is positive definite. 

It is easy to recognize that positivity of G is manifest once the adjoint scalars vanish 
<Pi = 0. Nevertheless, it is instructive and assuring if we consider the potential as well as the 
BPS equations keeping the adjoint scalars nonzero. 

To ascertain positivity of G we need to compute its eigenvalues. This is most easily done 
by looking at its determinant (We leave the derivation of this result to the Appendix IU0 : 



detG = 

Requiring det G > 0, we have 



m? m 2 , \ : 



\2 



2 ™2 



Ni f 1 \N 2 



P 2 \P 2 



(4-33) 



.^ + ^W 2 >0. (4-34) 
9192 V 91 92 ' 

In Appendix O we show that this condition is both necessary and sufficient to ensure posi- 
tivity of matrix G in Eq. ([4-240 . 

4.3. BPS equations 

Let us denote the codimension of the domain wall as y. Since we assume Lorentz in- 
variance for other dimensions, we obtain gauge field to vanish for component other than 

V- 
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The energy density 1-L for domain walls is given by 

H = l -G aP V y E a V y E p + \(G- l ) aP (r - c) a {r - c)p (4-35) 

2 

+ E Tr { (VyHiiPyHl + V y H iR V y H\ R + V y H} L V y H lL + V y H\ R V y H iR ) + V iF } , 
i=i 

where color-flavor indices a, (3 span all values as in Eq. (I4-4|) and we have incorporated color 
sector a = into the definition of matrix G for brevity. Accordingly, we have incorpo- 
rated the definition, (r — c)i i = r 1 *. Since there is no mixing of color sector with the rest, 
the inverse is calculated trivially and non-color part remains the same as in (14. ip . 

Now we observe that the mixing due to the cubic prepotential occurs only in the ki- 
netic term and potential of the vector multiplets. Moreover, they appear as G and G^ 1 
respectively. Therefore the cross term coming out of the Bogomol'nyi completion has no 
dependence on the metric G. This fact implies that the cancellation of cross terms to give 
topological charge goes through unaffected by the mixing of the vector multiplets. 

More explicitly, we obtain the Bogomol'nyi completion as 



H= l - [G^VyE^ + (r - c)}j (Cr 1 )^ \GpsVyS- 

2 

+ E Tr [fa 11 * + Zi R iL - H iL (#i + r^l 



Ni 

2 
|2 



+\V y H} L - SiSl + H\ L + mil Ni ) 
+ \V y H iR + SiH iR - H iR {$i - mil Ni ] 
+\V y H} R - Ei&h + H\ R (*< - ml Ni ) 

2 

- d y Tr [Ar< + ^ - TrniHlHiL - H\ R H lR - H lL H\ L + H iR H\ r ) 

i=l 

+ c a d y TrE a . (4-36) 

The last term gives the usual Bogomol'nyi bound and becomes the topological charge. The 
line before that is the total derivative which give vanishing contribution for an infinite line 
— oo < y < oo. 

BPS equations for iJ's and H's of hypermultiplets are 

V y H iL + EJIiL - Ha ($i + mil Ni ) = Q, (4-37) 

V y Hl - Ei&i + Hi (S, + ma*) = , (4-38) 

V y H iR + EiHtR - H iR (<2>* - rm l Ni ) = , (4-39) 

V y Hj R - E^l + H\ R (*« - ml Ni ) = . (4-40) 
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BPS equations for vector multiplets are 

G aP V v S p + (r-c) a = 0. (441) 

More explicitly, 

i 2 \ 

-±d y 2J > + Y Aj£ikmk $ A W y $ A i + r ' - c - = , (442) 
gf ^ mi m 2 



—V y E h + r h = , (443) 



1 

^ + y =V*a s ^ + s** = o . (444) 

e 2 ^ m im2 

We can easily solve the BPS equation for hypermultiplets, by using the moduli matrix 
approach. We define Si C , Si F and ipi as 

2i(y) + iW iy (y) = Sr\y)d y S ic {y) + \d y ^{y) , (445) 
$i(y) + iA iy {y) = S lF (y)d y S7 F \ y ) . (446) 

Then the hypermultiplets BPS equations (I4-37|) - (I4 40[) are solved by the constant moduli 
matrices Hf L and H® R 

H lL = e-WS^H^S^e™" , (447) 

H iR = e-^S-'H^e-^y , (448) 

where Si C , SiF £ SL(Ni,C). The hypermultiplet fields Hn and Hm do not contribute to 
domain wall solution and they are therefore vanishing. We write down ( I442p - (l444j) in terms 
of the gauge invariant fields 

f2 ic = S ic Sj c , f2 iF = Sj F S iF , Vi = l(^ + ^)- (4-49) 

The adjoint scalar fields of the vector multiplets are given by 

= \s; c \d y Q tc Q^) Stc + ld yVi , (4-50) 

^ = -\s^\d y n lF n^) S \ F . (4-5i) 

Also, we have 

DyS, = d y Ei + i [W iy , Et] = ^S^dyidyQiM^Sic + ^dfa , (*52) 

1 



Dy$i = d y $i + I [A iy , fy] = -^Sj F 1 dy(dyQ tF f2- F 1 )Sj F . (A^) 
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BPS equations for vector multiplets (I4-42|) - f!4-44|) can be now rewritten as the following 
master equations: 



-1\2 



i 2 + |^ Tr r A {Q i 

2gf y 2mim 2 L J y J 3 

= co t - e-^[(Hln^H^y + HhiftH&e-*™) Or 1 ] , (4-54) 

~2*-)y(dyf2i c fl ic ) 

= -e- rH ({H" lL n^Hf L e 2 ^y + H^H^^nr}) , (4-55) 



d y (d y (r] 1 m 2 - ■q2 l m 1 )d y n iF Q i p) 

m 1 m 2 

= -e~^((H^H^y + H&n^H&e-*"")^ 1 ) . (4-56) 
Here we have used a notation 

(X) = X - ^1n. (4-57) 

We make a comment about a possibility of additional moduli. At present we cannot say 
definitely if there are additional moduli other than the moduli matrices Hq's, since we cannot 
solve these master equations. We have several clues at hand. The BPS equations for domain 
walls and other solitons in gauge theories with scalar fields in the fundamental representations 
are in the Higgs phase where all the gauge symmetries are broken in the vacuum. In that 
situation, we learned that all the moduli are contained in the moduli matrix. On the other 
hand, instantons are solitons in the pure Yang-Mills theory without scalar fields, where 
gauge symmetry is unbroken in the vacuum. In this case, moduli reside in the BPS equation 
for gauge fields. In our present case, unbroken gauge symmetry SU(Ni) c+ v remains. This 
feature is indicative of additional moduli coming from the vector multiplet. 

Irrespective of the possible additional moduli, we can demonstrate that the BPS equations 
admit the coincident wall solution. Since the hypermultiplet parts are already solved as in 
(j4-47p - (|4-48p . our main task is to solve the master equations (!4-54p - (l4-56p associated to the 
vector multiplet. In order to solve them explicitly, we take strong gauge coupling limit 
gt,gi —> oo, where the master equations give just the algebraic constraints for f2 ic , Q iF and 
rji- In principle, they can be solved algebraically. Furthermore, Eq.( l4-34|) with the limit 
Qi — > oo tells us that positivity is maintained only if vanishes. In the following we will, 
therefore, consider a special point in the solution space where 



# f = 0, i = l,2 (4-58) 
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which implies from Eq. fl4-5ip that Qip are constant matrices. Then the differential equations 
fl4-55l) - fl4-56p reduce to the set of algebraic equations: 

H^H^y + H^H^e- 2m ^ = ^Q lc , (4-59) 

H^H^y + H^n-'H^e-^y = ^n iF . (4-eo) 

Notice, that for both sectors i — 1, 2 these equations are the same and do not couple to 
each other. We can, therefore, focus our discussion only on one sector, since all results are 
equivalent in both of them. So in the remaining discussion we will drop the index i from all 
fields. 

Now we consider moduli matrix for the coincident walls corresponding to the most sym- 
metric point of the moduli space. 

(H° L ,H° R ) = (1 N ,1 N ). (4-61) 

Eqs. fl4-59p and ( I4-60|) show that these two constant matrices commute and only the product 
can be determinec|*3 

&Q C Q F = — (e 2my + e- 2my ) 1 N . (4-62) 
c v 

Since we have chosen the matrices S c , Sp in SL(N, C), we find that det(f2 c f2 F ) = 1 and we 
can separate the U(l) part. 

e V = _ l e 2 my + e -2my\ = ^.gg) 

C v 

The U(l) part gives the usual domain wall solution. Without affecting the physical quan- 
tities, we can choose Q c = 1, S c = 1, and finally we obtain the coincident wall solution for 

dMU) with flMED 

= (27) = o , (4-64) 

n = log — ( e 2m &'-w) + e - 2m(y - yo) ) , (4-65) 
c 

S° = ^dyV — mtanh(2m(?/ — y )j , (4-66) 

f~c e m (y-y°h N 
Hl = \Ijj- / ( 4 - 67 ) 

o2m(y~y ) _|_ e -2m(y-y ) ) 



e m(j/-j/o)]_ 



H * = \ljf-, ~ ~ xl/2 - 



2m(y-y ) _|_ e -2m(y-y )^ 



*) This is due to the special choice of the moduli matrix in Eq. (l4-61|) . since Eqs. (14-471) and (|4-48[) imply 
that only the product SipSi c can enter into the physical fields such as Hi^^Hm. 
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Note that in this solution we restore a moduli parameter yo corresponding to the position 
of the coincident wall. A similar construction of domain wall solution works for the second 
sector (i — 2), besides the first sector (i = 1) given above. 

Let us note that the field-dependent gauge coupling function similarly to ( I3-33P is auto- 
matically obtained as a bosonic part of the Lagrangian specified by the cubic prepotential 
in Eq. (14-111) . Restoring the index i = 1, 2 for both of the domain wall sectors, and by using 
(14 -6 51) with (l4-50p . we finally conclude that the appropriate profile of the field-dependent 
gauge coupling function E Ql jm\ — E° 2 /m2, similarly to f!3-33j) is achieved. When we make 
(a part of) the global flavor symmetry as a local gauge symmetry, we can have several op- 
tions. Since the first flavor group SU(Ni) is generally different from the second flavor group 
SU(N2), we can naturally introduce two different gauge fields for the i = 1 and 2. This 
option leads to two decoupled sectors in the low-energy effective Lagrangian, which can only 
be coupled by higher derivative terms induced by massive modes. Another interesting option 
is to introduce a gauge field only for the diagonal subgroup of isomorphic subgroups of two 
different flavor groups, such as SU(N) G 3U(N X ), SU(N) G SU(N 2 ) with N < N u N 2 . This 
option is interesting in the sense that the massless gauge field exchange will communicate 
between two domain wall sectors. We hope to come back to these issues in near future. 

Let us make a few comments. First we have shown that the chiral model analyzed in 
section [3] can be extended to a supersymmetric gauge theory with eight supercharges and 
that the field-dependent gauge coupling function which is a clue for localization is naturally 
explained by taking the cubic prepotential. Second, there may be more moduli not contained 
in to (H®, H R ), which require further studies. Third, here we have presented a solution at a 
special point <& = 0. It would be interesting to consider the case or <P ^ 0, but in this case, 
we need to take a finite gauge coupling limit, on which we will investigate in future work. 

§5. Conclusions and discussion 

In this paper we have successfully localized both massless non-Abelian gauge fields and 
massless matter fields in non-trivial representation of the gauge group. We first considered 
a (4+l)-dimensional U(N) gauge theory with additional SU(N) L x SU(N) R x U(1)a flavor 
symmetry. We introduced the flavor gauge field for the diagonal flavor group SU(N) L+R , 
which is unbroken in the coincident wall background. The flavor gauge fields are localized on 
the wall by introducing the scalar-field-dependent gauge coupling function. Then we studied 
the low-energy effective Lagrangian and showed that massless localized matter fields interact 
minimally with localized SU(N)l +r gauge field as adjoint representations. Moreover, full 
nonlinear interaction between the moduli containing up to the second derivatives, was worked 
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out. The field- dependent gauge coupling function is naturally realized in supersymmetric 
gauge theories using the so-called prepotential. For this reason, we also explored bosonic 
part of M = 1 supersymmetric extension of our model. 

Main result of this paper is the effective Lagrangian (I3-49|) . The moduli field U appearing 
in the effective theory, is a chiral N x N matrix field like a pion, since it is a NG boson of 
spontaneously broken chiral symmetry. Other moduli in (13-491) . denoted by N xN Hermitian 
matrix x, has the physical meaning of positions of iV domain walls as its diagonal elements. 
We argued that the fluctuations of moduli field x, can develop VEV corresponding to splitting 
of walls, and the Higgs mechanism will occur as a result. Namely, the flavor gauge fields get 
masses by eating the non-Abelian cloud. Therefore, in this model, Higgs mechanism has a 
geometrical origin like low energy effective theories on D-branes in superstring theory. 

Amongst the possible future investigations, we would like to study non-coincident solution 
to further clarify this geometrical Higgs mechanism. 

We have noticed that our effective moduli fields resemble the pion in QCD. Similar 
attempts have been quite successful using D-branesP2> We believe that our methods can 
provide more insight in various aspects of low-energy hadron physics. We plan to explore 
this direction more fully in subsequent studies. 

In the discussion of supersymmetric extension of our model in section HI we employed a 
general setup where both sectors possessed their own domain wall solution, preserving the 
same half of the supercharges. But another alternative approach is also possible. We can 
consider a model, where different halves of supercharges are preserved at each sector (BPS 
and anti-BPS walls), and the SUSY is completely broken in the system as a whole. It has 
been proposed that the coexistence of BPS and anti-BPS walls gives the supersymmetry 
breaking in a controlled manner.^ 1 In our present case, BPS and anti-BPS sectors interact 
only weakly. If we choose flavor gauge field for each sector separately, we have only higher 
derivative interactions induced by massive modes. If we choose the diagonal subgroup of 
(subgroups of) each sector as flavor gauge group, we have a more interesting possibility of 
the massless gauge field as a messenger between two sectors. We plan to address this issue 
elsewhere. 

In order to construct a realistic brane-world scenario with the SM fields on the domain 
wall, we need the localization of fields in the fundamental representation of the gauge group. 
This is still an open problem and one of the priorities of our future investigations. In 
particular, the SM contains chiral fermions. Localization of chiral fermions is a particularly 
challenging problem. Anomaly associated with the chiral fermion is also an interesting issue 
to be addressed. We would also like to clarify these problems in subsequent studies. 

Two more issues remain to be addressed. First is the question of sign of gauge kinetic 
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term. In our present model, the positivity of the gauge coupling function is assured only 
when positions of walls are properly ordered (see Eq. (l3-480 ). namely only in a region of the 
moduli space, More economical models such as given in Ref JT2j) may not have such moduli 
and, therefore, the effective gauge coupling may be always positive. And lastly, as discussed 
in section HI we have not succeeded in exhausting all moduli in the supersymmetric extension 
of our model. We would also like to investigate these aspects in the future. 
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Appendix A 

Domain walls in the gauged massive CP 1 sigma model 

Here we consider the domain wall solutions in the gauged massive CP 1 sigma model. 
The model is obtained as the strong gauge coupling limit of a model similar to that we have 
studied in section I2T21 Namely, we start with the Lagrangian which has U(l) x U(l) gauge 
symmetry with two flavors 

C = -^Fmn? + ^(d M a) 2 - ^(Gmn) 2 + \V M H\ 2 - V, (A-l) 

V= j (\H\ 2 -v 2 ) 2 + \aH-HM\\ (A- 2) 

where H = (Hl,Hr). The covariant derivative is given by 

V M H = d M H + iw M H + ia M Hq, q = diag(g L , q R ). (A-3) 

The mass matrix is chosen M = diag(m, — m) as before. 

We next take the strong gauge coupling limit g — > oo of only one of the gauge coupling 
which results in the non-linear sigma model coupled to the other gauge field with the finite 
gauge coupling e. In the limit the gauge field % and the neutral scalar field a become 
Lagrange multipliers. After solving their equations of motion, we have 

w M = \V M HH\ a = \HMH\ (A-4) 
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where we have introduced the covariant derivative 

V M H = d M H + ia M Hq. (A-5) 

Plugging these into the original Lagrangian at g — > oo, we get the gauged massive CP 1 
sigma model 

eg-** = -^(Gmn) 2 + V M HPV M H^ - HMPMH\ (A-6) 



with the projection operator 



P = 1 - -riZtff. (A-7) 



V 2 



As before, let us rewrite this Lagrangian with respect to the inhomogeneous coordinate 

H= 7 =JL== (1,0), (f) E C. (A-8) 
V 1 + 101 

Then the charge matrix should be chosen as 

g = diag(0,l), (A-9) 

which leads to a natural expression that the complex scalar field <fi has the U(l) charge 1 for 
the gauge field oj^: 

V M H = - 2(1+ [ 0|2)3/2 (^|0| 2 , 09 M |0| 2 - 2(1 + |0| 2 )P M 0) , (A-10) 
V M (l>=(dM+iaM)<P- (A- 11) 



Plugging these into Eq.( lA-6p . we finally get the Lagrangian 

1 2 2 |£>m</>| 2 2 4m 2 |0| 2 

£ S -*oc = -T-^(W) + U JZ , i ,, 2 , 2 - U 73 ; . ,. 2 , 2 . (A-12) 



Let us next consider a domain wall solution in this model. We assume all the fields 
depend on only the extra-dimensional coordinate y. Then the four dimensional components 
of the Maxwell equation 



,V M (f)(j)* - <f)V M - 

:i + i0i 2 ) 2 ' 



o N g NM = z e v - ,j~ 2 y , (a-is) 

can be immediately solved by 

a M = 0, (i = 0,l,2,3. (A-14) 
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The fifth component is 



= ie v 



2„.2 



V y 66* 



1 + 



12^2 



(A-15) 



Now the Hamiltonian reduces to the following form 

,2 



H 



V 



;i + |0|2)2i' y 

,2 



V v 6r + 4m 



2 1 / 1 2 



V 



> 2mv 2 



d i 



Vyd) + 2mc 



-2md y \(j>f) 



(A-16) 



Thus the reduced Hamiltonian is minimized when the following first order equation is satisfied 

V y( j) = -2m<f>. (A-17) 
Since the mass parameter m is real, Eq. (1A-15j) is also satisfied. Let us take the gauge where 



a y = 0. 



Then we have the explicit domain wall solution 

i _ Q2 e ~2my q2 _ e 2ia+2my 



(A-18) 



(A-19) 



This is completely the same as the domain wall solution given in Eq. ( I2-28P in the ungauged 
massive CP 1 sigma model. 

The final step is to obtain a low energy effective theory on the domain wall. The effective 
Lagrangian is given by 



£ 



eff 

S->-oo 



dy 



dy 



1 ^2 I 



,2 2 



4e2 (W 2e 2^w) +V (1 + |0|2)2 

v 2 (m 2 (d,y,) 2 + (V,a) 2 



1 - ^2 1 (n ^2 



(A-20) 



where we have promoted the moduli parameter yo, a to the fields 2/0(2^), o;(a; /1 ) on the wall, 
and we have introduced the covariant derivative 



(A-21) 



where a is the function of the (3+l)-dimensional coordinate x^. Assuming a M to be y- 
independent (zero mode), we finally obtain 



'eff 

-g->-oo 



1 



m 



-4-2^) + - (m^yoY + (X» 



(A-22) 
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Thus we find that the gauge field a^(x) absorbs the scalar field a(x) to become massive 
via the Higgs mechanism. Since that the U(l) gauge field a M is massive in the effective 
Lagrangian, we have to integrate it out according to the spirit of the low energy effective 
theory. 



Appendix B 

Effective Lagrangian on the domain wall 



In this appendix we derive our main result (13-41 p of the effective Lagrangian for the 
gauged Chiral model introduced in §3J 

B.l. Compact form of gauged nonlinear model 

Starting from the Lagrangian using the Einstein summation convention for a = {L, R} 



off 



dyTi 



D^H a D»Hi - v 2 WAV» 



with the constraint 



H a H\ = v 2 l 



N 



(B-l) 



(B-2) 



we first eliminate the gauge fields to obtain a simple expression for gauged nonlinear 
sigma model. Gauge fields are given by equations of motion as 



2v 2 



D,H a H\ - H a D,Hl 



and 



(B-3) 



(B-4) 



The effective Lagrangian (IB- ip should also contain kinetic term for gauge field A^, but we 
will not explicitly write it here, for brevity. Eq. (IB-ll) can be further simplified by using the 
following identities 



H a D^Hl = d,(H a Hl) - D,H a Hl , 
Hlt)^H b = -D,H a H b + V^H ab , 



where 



T^pMab = d^Hah + i [Ap, H ab ] 
After some algebra we find: 



H a b = HlH h 



(B-5) 
(B-6) 

(B-7) 
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Plugging above expression back into the (IB-lj) we arrive at: 



-y poo 
£pff = :r^r 



2t> 2 



dy Tr 



'D IJr H ab T> tJ 'H ba 



(B-8) 



B.2. Effective Lagrangian 

Now we are ready to compute effective Lagrangian. Using a solution (with y = mylN—x): 



H=(H L ,H R ) 
Our new fields H a t> are given as: 

Hll 
Hlr 
Hrr 



oV/2 



V2 Vcosh(y) ' y/2 y/cosh(y)J 



(B-9) 



2 cosh(y) ' 



2 cosh(y) 
v\ v _ e- 



RL ' 



2 cosh(y) 
It can be checked, that (1B-8|) is given as: 

2 /-oo r p?/ 



(B-10) 
(B-ll) 
(B-12) 



cff 



4 



P, 



1 1 



+£/ f p M ^ 



— oo 



(_ cosh(y) cosh(^) cosh(y) cosh(y) 



cosh(y) 



iflw (u 



cosh(y) 



cosh(y) ' cosh(y) 



cosh 



(B-13) 



In the following we would like to carry out the integration over the extra-dimensional coor- 
dinate y. This can be done in two steps. First, we must factorize all quantities depending on 
y (or on y) to one term inside the trace, effectively reducing our problem to fit the following 
form: 



dy Tr f{mylN — x)M 



(B-14) 



where M is some matrix, independent of y and / some function. In the second step we 
diagonalize x: 

x = P -1 diag(Ai,...,AAr)P, 
and use the fact that f(P~ 1 yP) = P^ 1 f(y)P. This transformation leads to 



dy Tr 



f(myl N - diag(A l ))PMP- 1 j = / dyJ2f( m V ~ ^(PMP" 1 



8=1 
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For every term in the sum we can perform substitution y = my — Aj. The key observation 
is that in each term the integration will be the same and independent on a particular value 
of Aj. Thus we arrive at an identity 

1 



dyTr f(y)M 



rn 



Tr(M) / dyf(y). 



(B-15) 



It appears as if we just made a substitution y = yljv- This is possible, of course, only thanks 
to the diagonalization trick and properties of the trace. In the subsequent subsections, 
however, we will refer to this procedure as if it is just a 'substitution', for brevity. Let us 
decompose the effective Lagrangian flB-13j) into three pieces 



^eff Tx ~\~ Til Tmixed 

and see the outlined procedure for each term. 



(B-16) 



B.2.1. Kinetic term for U 

First, let us concentrate only on terms containing double derivatives of U, which we 
denote 7u' 



% = V - dyTrlV^WWU 



1 



cosh 



2 +V^U 



_cosh(y) 



cos 



m } 



where we have used the fact that inside commutator it is possible to freely interchange 



cosh(^) cosh(y) 

since the difference is just a constant matrix. In this way we made 7u manifestly invariant 
under exchange y — >■ —y. 

Since in the first factor of Tu all y-dependent quantities are on the right side, we can, 
according to our previous discussion, make use of the identity (IBT5j) and carry out the 
integration: 



— / dyTi 

^ -' — oo 



V u^v^u - 

cosh 2 (y) 



Tr 

2m 



For the second term, however, we first use the identity: 

[f(A),B} = Y,^ k A (B)f k \A) 



k=l 



where Ca(B) = [A,B] is a Lie derivative with respect to A. Thus 



cosh(y) 



E 

k=i 



ci(u^v»u) 



\ (*) 



cosh(y) 



(B-17) 
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Now all ^-dependent factors are standing on the right and we can formally exchange y — >■ y. 
The summation can be carried out to get: 



°° (l\k / y \ (fc) 



cosh(y) 

The formula for Tu now reads: 



^ 4m 



dy Tr 



cosh(y - £ 4 ) 



— r • 

3 *—a: 



(U*WU) - 



cosh(y) 



U*WU. (B-18) 



cosh(y - C x ) cosh(y) 



(U^U)V^U 



(B-19) 



Since we started with 7[/ invariant under the transformation y — > —y, we should take only 
even part of the above formula (under exchange C £ — > as the final result: 



^ U Am 



dy Tr 



cosh(£i 



cosh(y - C x ) cosh(y) 



-(U^U)V^U 



(B-20) 



Now we can carry out the integration using primitive function 

dy 1 . 1 



In 



cosh(|/ — a) cosh(y) sinh(o;) 1 — tanh(a) tanh(y) 

Therefore we obtain the result to all orders in x as: 

1 , /l + tanh(£ £ ) 



y ^ 

4m 



T^fjU ^ ^^[Cx) ^ \l — tanh(£ £ ) 



(B-21) 



Performing the Taylor-expansion of the function 



1 



In 



1 + tanh(x) 



2x & 



Ax 



10 



2.x 2 2x 4 Ax 6 

tanh(x) " V 1 - tanh(x) y + IT ~ lo" + 945 ~ 4725 4 93555 



+ 0(a; 



12\ 



(B-22) 



we can easily read off coefficients of terms beyond the leading one. For example, the first 
three terms reads: 



% ~ 2m " V" " 



v 

6m 



- ^-Tr ( [x, [x, U^V.U] ] [x, [x, Wrfu] ])+... (B-23) 



B.2.2. Mixed term 

Mixed term between x and U is given by 

„,2 roo 



mixed 



dyTrtlJ^U 



1 



cosh(y) ' cosh(y) 



cosh(y) ' cosh(y) 
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With use of the identity (1B-17j) and 



X 



k=0 



f( k+l \x) 

' (k + iy. 



(B-24) 



one can prove the following: 



[f(x),V,g(x)) =J2ri C t\V,x)[(mg(x) 



n=2 



(«) 



_ f( n ) 



fW(x)g(x)-f(x)g( n \x) 



We can use this result to factorize all y-dependent quantities to the right and make the 
substitution y = yl^\ 



,2 r°° 



mixed 



4m I d ^ 2 



Tr 



n=2 



IV. 



(n) 



X 



1 



(fl) 



cosh(y) / cosh(y) \ v cosh(y) / cosh(y) \cosh(y) J cosh(y) 
Now we are free to perform summation and integration to obtain: 

cosh(£ £ 



Tmixed „ Tr 

2m 



L x smh(>L; 



Performing the Taylor-expansion of the function 



cosh(x) — 1 f 1 + tanh(x) 
In 



oc oc 



x — — — -)- 



1 — tanh(£^ 



17x 7 31s 9 



(B-25) 



+ 



+ 0(x Li ), (B-26) 



xsinh(x) \1 - ta,nh(x) J 12 ' 120 20160 ' 362880 

we can easily read off coefficients of terms beyond the leading order in the series expansion: 



Tfjiixed r. Tr 

2m 



24m 



-Tr 



WV^U [x, [x, [x, V^x] 



+ 



(B-27) 



B.2.3. Kinetic term for x 

Kinetic term for x is given by 

,2 poo ( | 



% = - 

I X 



-2> 



1 



4 J_ OQ ° \ M cosh(y) cosh(y) M cosh(y) cosh(y) J 
We are going to need the identity 

. i-iy 



(B-28) 



VJ(x)V»g(x)\ =J2 y -^ L TrUr\V tl x)V»x 



n=2 



(n) 



f(x)g(x)) -f^(x)g(x)-f(x)g^(x) . (B-29) 
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With the aid of this we arrive at 



v 

4m 



00 n=2 
J) 



X 



(n) 



737 + 



(n) 



1 



1 



cosh(?/) / cosh(?/) \cosh(?/) / cosh(?/) \cosh(?/) / cosh(y) 

where we again employed diagonalization trick and identity (IB- 15[) . Let us carry out the 
summation and the integration to obtain: 



% = Tr 

2m 

leading to the power series: 



In 



1 + tanh(£.; 



v £ cosh(£ a ) 
M £|sinh(£ :? .) "\l — tanh(£ £ ) 



% = Tr 

2m 



V^xV^x 


v 2 m 


+ Tr 




24m 



+ 



(B-30) 



(B-31) 



Putting all pieces together as £ e ff = % + 7u + l~ m i xe d, we obtain our final result (13-411) . 



Appendix C 

Determinant of G 



In order to calculate determinant of matrix G (I4-24|) we will use the following recurrence 
formula, which relates determinant of a symmetric matrix AA of rank N + 1 to determinant 
of its iV x N submatrix M: 



a u 
u M / 

(a-v^M^u) detM. 



M 

detM 

After double application of formula (\C-2\i we get 

detG = 

/ \ V \ 



1 



0.1 \ ' mi ' mi / 



V 



/ 



A^Ai 
mi 



ni2 



\ 



-1 



Go 



_AL(^4i 



UJVl / 1 \ ^2 

xK x 



(C-l) 
(C-2) 



(C-3) 
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where 



Gi 



Go 




tS a 2 b 2 



m 2 





(C-4) 



(C-5) 



Inverse of G\ is given as 



/ 



ft 



cr 1 



— x 



1 - %m\<f> 2 



e 2 nng2 



2$A 2 



\ 



eimigl$ Bl e\5 MBl - m\e\g\$ MBx 
\e 2 m 1 gl<P B2 mje 1 e 2 gl$ B2 $ Al 



in 



leie 2 gl$ Bl $ M 



e 2 5 A * B * - m 2 e 2 g 2 S A ^ J 



where we have used (I4-26M-29]) . 

Straightforward calculation leads us to 



detG 



ml<P 2 



1 - glm\^> 2 - 



- m\<P 2 



^2 



1 x Ni , 1 \ N 2 



After multiplying both brackets we obtain the result (I4-33P : 



(C-6) 



(C-7) 



detG 



2 - » 

91 92 



A'l 



i / 1 \ N 2 



(C-8) 



Next we would like to find condition, which ensures positive definiteness of G. In other 
words, we require that all eigenvalues of G are non-negative. We can easily turn (I4-33P into 
characteristic equation by replacing g^ 2 , g^ 2 , ej~ 2 , 2 with g^ 2 — A, g^ 2 — A, e\ 2 — A, 2 — A. 
However, since §> 2 consists of terms proportional to either e\ or e 2 , instead of ej~ 2 or e^~ 2 , 
we should first multiply the term in the square bracket by a factor e^e^ 2 ■ Then, after the 
replacement and denoting <P 2 = <& Ai (I? A \i = 1,2, we obtain a characteristic equation of the 
forth order: 



A 



L v #r 



1 

¥2 



A 



times the factor 



1 



-A 



- - 



1 



A 



91 



+ 



9l 



A(m 2 + m\) 



$2 



A 



2 2 
m\m2 



A 



iVi-l 



2V 2 -i 



(C-9) 



(C-10) 
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which clearly leads to positive eigenvalues. Expanding the brackets we obtain explicit coef- 
ficients of the characteristic polynomial: 

( 1 1 1 1 1 1 m l+ m L\2^2 \2^2^\2 

+ 1 T9T9 + To 9 + „o 9 + „ 9 9 + „ 9 9 H 9- 9~ 9 5 (Ai^i + ^2^2) ) A 

V #l#2 #l e l #l e 2 #2 e l #2 e 2 e l e 2 K m 2 7 

^ + ^ + 61 + 62 _ ~ 2 A 2 <£ 2 + Ap> 2 _ m|+m| ~ 2 x ^ + 1 - g 2 $ 2 
g\gle\e\ m\m\g\g\ e\e\ J g\gle\e\ ' 

In order to see the non-negativeness of eigenvalues it is not necessary to solve the char- 
acteristic equation. Generally speaking, characteristic equation of a real symmetric matrix 
can be always put into the form 

(A-A 1 )(A-A 2 )...(A-A JV )=0, (C-12) 

where all roots Ax,... Xn axe real numbers. Multiplying all parentheses we see that the 
coefficients c& of characteristic polynomial are given by the sum of all possible A;-tuples of 
A's with alternating sign: 

(A-A 1 )(A-A 2 )...(A-A JV ) 

= X N -(£ K)^- 1 + QT A,A,) A^ 2 - . . . + (-1)*% . . . Ajv 

i i>j 

N 

= ^2{-l) i c i X N - Ji . (C-13) 

i=0 

The positivity of all coefficients Ck turns out to be equivalent to the positivity of all eigen- 
values Afe. To ensure positivity of the eigenvalues, we now demand that all terms in (ICTip 
inside brackets are positive. This gives us three conditions: 

-2^2 , -2 2 , -2 2 , -2 2 , -2 2 , 2 2 
9l9 2 + #l e l + #l e 2 + 92 e l + 92 e 2 + e l e 2 

-g\gl{m\ + m 2 )(e 2 <£ 2 + e\$\) > , (C-14) 

gl + 9l + el + e 2 2 - ~g\e\$l + e\$\) - g\g\{m\ + m\)& > , (C-15) 

l-f$ 2 >0. (C-16) 

These can be put into the convenient form: 

l-f 2l $l-g 2 22 <Pl>0, (C-17) 
l-g 2 n S 2 -g 2 12 3 2 2 >0, (C-18) 
l-g 2 3l-g 2 3 2 2 >0, (C-19) 
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where 

~ 2 = glgl(jn\ + m 2 2 )e 2 2 

9i9l + 9 2 A + g\e\ + g\e\ + g 2 e 2 + e 2 e 2 ' 

~ 2 = gjgjjmj + m 2 2 )ej 

fe 2 + 9 2 i4 + g\e\ + g\e\ + g\e\ + efe| ' 

~ 2 _ <7 2 e| + g?flf(m? + mp 

" n " gf + gi + ej + el ' (C ' 22) 

~ 2 _ g 2 el + glgl(m\ + ml) 

9x2 ~ gl + gl + el + el ■ {L ' 26) 

We are going to argue, that the last condition 1— g 2< P 2 > is the strongest one and, therefore, 
only one important. This can be true if and only if parameter g 2 is always greater then gf l: 
gf 2 , g\\ and g 2 2 f° r an possible values of involved parameters. This is indeed so. Let us 
demonstrate this fact by showing, for example 

f > ~g 2 n ■ (C-24) 
Multiplying both sides by g\ + g\ + + e 2 and expanding our notation, we get 

[m\g\ + m\gl){g\ + g 2 2 + e\ + e\) > e\{m\g 2 2 + m 2 2 g\) + g\g 2 2 {m\ + m 2 2 ), (C-25) 
leading to 

™ 2 ^ 2 2 + e 2 ) + mlg 2 x {g{ + e?) > . (C-26) 

Last line is obviously always true. In the same way, one can show that g 2 > gf 2 , g 2 > g 2l 
and g 2 > ^| 2 - This proves our claim, that condition f!4-34j) is both necessary and sufficient 
to ensure positivity of the matrix G. 
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